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Abstract. - Atoms, propagating across a detuned standing laser wave, can be scattered in a 
chaotic way even in the absence of spontaneous emission and any modulation of the laser field. 
Spontaneous emission masks the effect in some degree, but the Monte Carlo simulation shows that 
it can be observed in real experiments by the absorption imaging method or depositing atoms on 
a substrate. The effect of chaotic scattering is explained by a specific behavior of the dipole 
moments of atoms crossing the field nodes and is shown to depend strongly on the value of the 
atom-laser detuning. 



• Introduction. — The deflection of an atomic beam 
(— | at a laser standing wave (SW) is explained by the dipole 
O forces which are well described by the classical atom-field 
L - — interaction model [TJ|2]. The ability of a SW to diffract, 
j focus or splitting an atomic beam [3] has been used for 
a variety of applications including atom microscopy, in- 
■^j" terferometry, isotope separation, and optical lithography 
C ■ @HS]- On the other hand, cold atoms are ideal candidates 
to test fundamental principles of quantum physics includ- 
ing the phenomenon of dynamical chaos at the microscopic 
level that is known as a kind of random-like motion in 



a deterministic environment [7T416j . Dynamical chaos is 



characterized by exponential sensitivity of trajectories in 
• »the phase space to small variations in initial conditions 
^ and / or control parameters. 

It has been proposed in Ref. [7] to study quantum chaos 



X 

j^^nd the corresponding effect of dynamical localization 
placing cold atoms in a far-detuned SW with a periodic 
kick- like modulation of positions of the SW nodes. A num- 
ber of experiments [9] have been carried out in accordance 
with this proposal. At large detunings, the atoms are not 
excited being quantum analogues of classical kick rotors. 
Since those experiments on atom optics realization of the 
5-kicked quantum rotor, cold atoms provide new grounds 
for experiments and theory on quantum chaos. 

It has been shown in Ref. [TU] that even a single-pulse 
far-detuned SW can induce chaos in atomic motion. For 
sufficiently large detuning, the excited state amplitude can 



be adiabatically eliminated [7], leading to a Hamiltonian 
with an external degree of freedom only. The correspond- 
ing equations of motion for an externally modulated non- 
linear pendulum constitute the well-known model with one 
and half degree of freedom that can be chaotic under some 
conditions. The other possibility is to induce chaos in spin 
degree of freedom of atoms periodically kicked by applying 
short magnetic field pulses [TT]. That is the one and half 
degree of freedom model of a kicked top. 

In difference from those and other papers on the re- 
lated topic, we consider the physical situation with com- 
paratively small detunings and should take into account 
a coupling between external and internal atomic dynam- 
ics, leading to a model with three degrees of freedom. It 
will be shown in the present paper that in this case one 
needs no modulation or any other perturbation of the SW 
to induce chaotic internal and external dynamics of atoms 
crossing the SW laser field. 

Near the atom-field resonance, when the interaction be- 
tween the internal and external atomic degrees of free- 
dom is intense, there is a possibility to create conditions 
for chaotic behavior without any kicking and modulation 
[15UT7] . If so, it is open the way to test the novel regime 
of atomic motion caused by the peculiarities of the dipole 
force in the strong coupling regime. In the semiclassi- 
cal approximation, atom with quantized internal dynam- 
ics is treated as a point-like particle with the rlamilton- 
Schrodingcr equations of motion constituting a nonlinear 
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dynamical system [T5HT5] . In a certain range of the atom- 
field detunings, a set of atomic trajectories becomes expo- 
nentially sensitive to small variations in initial quantum 
internal and classical external states or/and in the con- 
trol parameters. Hamiltonian evolution is a smooth pro- 
cess that is well described in a scmiclassical approxima- 
tion by the Hamilton- Schrodinger equations. The prob- 
lem becomes much more complicated because of sponta- 
neous emission of atoms with a specific shot quantum noise 
acting in a dynamical system which is chaotic in the ab- 
sence of noise. A number of nonlinear Hamiltonian and 
dissipative effects have been found numerically and ana- 
lytically near the resonance including chaotic Rabi oscil- 
lations, chaotic atomic transport, dynamical fractals, and 
Levy flights Q3H2]. 

The main aim of the paper is to demonstrate theoreti- 
cally and numerically that the new type of atomic diffrac- 
tion at a rigid SW without any modulation, chaotic atomic 
scattering, can be observed in a real experiment. The 
scheme of such an experiment is shown in Fig. [1] with a 
beam of atoms crossing a SW laser field. One either mea- 
sures a spatial atomic distribution after the interaction by 
the absorption imaging technique or measures an atomic 
distribution on a silicon substrate in the far-field zone. 
The results are expected to be different depending on the 
value of the atom-field detuning. The distribution is ex- 
pected to be comparatively narrow at those values of the 
detuning at which atomic scattering is regular (r.s. distri- 
bution in Fig. [T]) or wide at the detuning values providing 
chaotic scattering of atoms due to their chaotic walking 
along the SW (c.s. distribution in Fig.[T]). 

atomic beams 



the peak value. The longitudinal velocity of atoms, v z , 
is much larger than their transversal velocity v x and is 
supposed to be constant. Thus, the spatial laser profile 
may be replaced by the temporal one. The Hamiltonian 
of an atom in the ID SW field can be written in the frame 
rotating with the laser frequency w/ as follows: 



M7 exp[-(t - -cr t ) 2 /of] (<7_ + d + )coskfX —a+a^, 

(1) 

where a± >z are the Pauli operators for the internal atomic 
degrees of freedom, X and P are the classical atomic po- 
sition and momentum, T, us a , and flo are the decay rate, 
the atomic transition and maximal Rabi frequencies, re- 
spectively, and ot = f/v z , i.e., 2>o~ t is the transit time. 
The wavefunction for the electronic degree of freedom is 
= o(*)|2) + 6(t)|l), where a = A + ia and b = B + i/3 
are the probability amplitudes to find the atom in the ex- 
cited, |2), and ground, 1), states, respectively. 





Fig. 1: Schematic representation of the proposed experiment on 
regular (r.s.) and chaotic (c.s.) atomic scattering at a Gaussian 
standing laser wave. 

Results. 

Main equations and the regimes of atomic motion. 
A beam of two-level atoms in the z direction crosses 
a SW laser field with optical axis in the x direction 
(Fig.[l}. The laser field amplitude has the Gaussian profile 
exp[— (z — z n ) 2 /r 2 ] with r being the e~ 2 radius at the laser 
beam waist. The characteristic length of the atom- field in- 
teraction is supposed to be 3r because the light intensity 
at z = zq ± 1.5r is two orders of magnitude smaller than 



Fig. 2: Finite-time Lyapunov exponent A vs atom-field detun- 
ing A (in units of the Rabi frequency $lo) and initial atomic 
transversal momentum po (in units of the photon momentum 
Kkf) at the normalized recoil frequency ui r — 1CP 3 and 7 = 0. 

In the semiclassical approximation, atom with quan- 
tized internal dynamics is treated as a point-like particle 
(the transversal atomic momentum p is supposed to be, 
in average, much larger than the photon one, hkf) with 
the equations of motion written for the real and imaginary 
parts of the probability amplitudes 

3 

x = ui r p, p = — 2cxp[— (r — -a T ) 2 / a 2 .](AB + aft) sinx, 

11 3 

A = -{uj r p 2 - A)a- - 7 A-exp[-(T- -a T ) 2 /a 2 } f3 cos x, 

113 

a = --(uj r p 2 - A) A - -7a + exp[-(r- -a T ) 2 /a 2 ]B cosx, 

1 3 
B = -(uj r p 2 + A)P - exp[-(r - -a T ) 2 /a 2 ]a cosx, 

$ = -\{u r p 2 + A)S + exp[-(r- '^ct t ) 2 / a 2 } A cosx, 

(2) 
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Fig. 3: Evolution (r is in units of J7q 1 ) of the atomic dipole- 
moment component u = 2(AB + a/3) in (a) the chaotic (A = 
0.2) and (b) regular (A = 1) regimes of atomic motion. 



Fig. 4: Trajectories in the real space for 50 atoms without 
spontaneous emission. Hamiltonian (a) chaotic (A = 0.2) and 
(b) regular scattering (A = 1). The atomic position x is in 
units of the optical wavelength. 



where x = kfX and p = P/hkf are atomic center-of-mass 
position and transversal momentum, respectively and dot 
denotes differentiation with respect to the dimensionless 
time t = £lot. The recoil frequency, uj t = fifc^/m Q Oo -C 1, 
the atom-laser detuning, A = (cof — uj a )/£la, the decay 
rate 7 = T/Qq, and the characteristic interaction time, 
<j T = rQ,Q/v Zl arc the control parameters. 

Without the decay, 7 = 0, Eqs. @ constitute the non- 
linear Hamiltonian dynamical system with three degrees of 
freedom describing an atom moving in the six-dimensional 
phase space. The simple way to know how complicated 
this motion may be is to compute the quantitative mea- 
sure of chaos, maximal Lyapunov exponent characteriz- 
ing the mean rate of the exponential divergence of ini- 
tially close trajectories. Because of a transient character 
of chaos we compute the finite-time Lyapunov exponent A, 
i.e. the value of the exponent at the moment when atoms 
leave the interaction region. The result of computation 
with Eqs. © at the given value of the recoil frequency, 
uj r = 10~ 3 , and zero decay rate in dependence on the de- 
tuning A and the initial atomic transversal momentum po 
is shown in Fig. [2] Color codes the magnitude of the finite- 
time Lyapunov exponent. In white regions the values of 
A are almost zero, and the internal and translational mo- 
tion is regular in the corresponding ranges of A and po- 
In shadowed regions positive values of A imply unstable 



motion. 

The scheme of the proposed experiment in Fig. Q] re- 
sembles the scattering problem with particles entering an 
interaction region along completely regular trajectories 
and leaving it along asymptotically regular trajectories 
[23 28] . However, in difference from the standard chaotic 
scattering with a nonattractive fractal invariant set exist- 
ing over an infinite time, this process may be interpreted 
as a transient chaos or a finite-time chaotic scattering. 

There are three possible chaotic regimes of the center-of- 
mass motion along the SW optical axis [17l[T8] . In depen- 
dence on the initial conditions and the values of the control 
parameters, atoms may oscillate chaotically in wells of the 
optical potential or move ballistically over its hills with 
chaotic variations of their velocity. Chaotic motion be- 
comes possible in a narrow range of the detuning values, 
< |A| < 1. AtA = 0, the synchronized clectric-dipolc 
component, u = 2(AB + aj3) becomes a constant. That 
implies the additional integral of motion in the Hamilto- 
nian version of Eqs.([2]) and the regular motion with A = 0. 
Far from the resonance, at |A| > 1, the motion is again 
regular both in the trapping and flight modes. That spec- 
ulation is confirmed by the Lyapunov map in Fig. [5J 

Moreover, there is a specific type of motion, chaotic 
walking in a deterministic optical potential, when atoms 



p-3 



S. V. Prants 




1000 T 1500 



Fig. 5: The same as Fig. 2] but in the momentum space. 



can change the direction of motion alternating between 
flying through the SW and being trapped in its poten- 
tial wells. We would like to stress that the local insta- 
bility produces chaotic center-of-mass motion in a rigid 
SW without any modulation of its parameters in differ- 
ence from the case with periodically kicked and far de- 
tuned optical lattices [71I5HTT]. The trivial time depen- 
dence in the Hamiltonian (|2) cannot produce chaotic mo- 
tion, it simply accounts for a modulation of the interac- 
tion of atoms with a Gaussian laser beam. Even if the 
atoms would cross an absolutely homogeneous (in the z- 
direction) laser beam there would be under appropriate 
conditions chaotic atomic center-of-mass motion in the 
transversal x-direction. 

Chaotic walking occurs due to the specific behavior of 
the Bloch-vector component, u, of a moving atom whose 
shallow oscillations between the SW nodes are interrupted 
by sudden jumps with different amplitudes while atom 
crosses each node [TS]. We illustrate in Fig. [3] the behavior 
of the u component with different values of the detuning 
A = 0.2 and A = f at which the atomic motion in ac- 
cordance with the A-map in Fig. [5] is chaotic and regular, 
respectively. The time of the atomic interaction with the 
SW field is estimated to be 3a T = 1200. So, the jumps of 
the u variable (if any) disappear after that time in Fig. [3] 
It follows from the second equation in the set © that 
jumps in the variable u = 2(AB + a/3) result in jumps of 
the atomic momentum while crossing a node of the SW. If 



the value of the atomic energy is close to a separatrix one, 
the atom after the corresponding jump-like change in p can 
either overcome the potential barrier and leave a potential 
well or it will be trapped by the well, or it will move as 
before. The jump-like behavior of u is the ultimate reason 
of chaotic atomic walking along a deterministic SW. 

It is easy to estimate the range of initial momenta at 
which atoms are expected to change their direction of mo- 
tion or move ballistically. At small detunings, |A| <C 1, the 
total energy consists of the kinetic one, K = ui r p 2 /2, and 
the potential one, U = it cos a;. If K(t = 0) > J7 max | = 1, 
then the atom will move ballistically. This occurs if the 
initial atomic momentum, satisfies to the condition 
Po > y/2/uj r > 44. If the initial conditions are chosen 
to give < K(t = 0) + U(r = 0) < 1, the atoms with 
< po < 44 are expected to perform a chaotic walking at 
positive A. 
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Fig. 6: Scattering of 50 spontaneously emitting atoms at the 
SW with the decay rate 7 = 0.05 and the same other conditions 
as in Fig. H (a) Chaotic (A = 0.2) and (b) regular (A = 1) 
regimes. 

Hamiltonian chaotic scattering. Let us consider a spa- 
tially uniform and previously focused beam of atoms cross- 
ing a Gaussiam laser beam. The position and momen- 
tum distributions of atoms are measured after interaction 
with the SW field. We predict that those distributions 
would be much broader at those values of A at which 
one expects chaotic walking to occur. Firstly, we per- 
form simulation with a negligible probability of sponta- 
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Fig. 7: The distributions of 10 4 lithium atoms (a) without 
and (b) with spontaneous emission at r = 1000 under the 
conditions of chaotic scattering at A = 0.2 (bold curves) and 
regular scattering at A = 1 (dashed curves). 



ncous emission. To be concrete let us take calcium atoms 
with the working intercombination transition 4} So — 4 3 Pi 
at A a = 657.5 nm, the recoil frequency i> rcc ~ 10 KHz, 
and the lifetime of the excited state T sp = 0.4 ms. Taking 
the maximal Rabi frequency to be £Iq/2tt = 2 • 10 7 Hz, 
the radius of the laser beam r = 0.3 cm, and the mean 
longitudinal velocity v z = 10 3 m/s, the interaction time is 
estimated to be 0.9 ms. The normalized recoil frequency 
is u! r = 47r^ rec /r2o = 10 -3 and oy — 400. 

We numerically solve the equations of motion ([2]) with 
7 = at two values of A corresponding to chaotic and 
regular regimes of the center-of-mass motion. In accor- 
dance with the Lyapunov map in Fig. [3J behavior of the 
Bloch component u in Fig. [31 and the simple estimates 
given above, we expect the chaotic scattering of atoms at 
A = 0.2 and their regular motion at A = 1. Trajectories 
in the real and momentum spaces for 50 atoms with the 
same initial momentum, pq = 10, and initial positions in 
the range — 7r/10 < x < w/10 are shown in Figs. [Hand 
5, respectively, at the fixed value of the recoil frequency 
Lu r = 10~ 3 . The upper panels in Figs. [4] and 5 illustrate 
the broad distributions of the atoms in the x and p spaces 
in the regime of chaotic scattering that contrasts strictly 
with those obtained in the case of the regular scattering 
at A = 1 (Figs. HJd and 5b). In order to create narrow 



atomic beams, one may use a pair of light masks. The 
first SW with a red large detuning (A < 0) splits the ini- 
tial atomic beam into a number of narrow beams with the 
widths much smaller than the optical wavelength which 
then cross the second SW. The method for creating nar- 
row wave packets in the nonadiabatic regime of scattering 
has been proposed in Ref. [5§] ■ 

Dissipative chaotic scattering and simulation of a real 
experiment. We have illustrated in Figs. [4] and 5 the 
Hamiltonian chaotic scattering that may occur in the ab- 
sence of any losses. To simulate trajectories of sponta- 
neously emitting atoms we use the standard stochastic 
wave- function technique (see, for example, [3"0"H3"2"] ) for 
solving Eqs. @. The integration time is divided into 
a large number of small time intervals St. At the end 
of the first one r = n the probability of spontaneous 
emission, s\ = 7<5r|a Tl | 2 /(|a Tl 2 + \b Tl \ 2 ), is computed 
and compared with a random number e from the inter- 
val [0,1]. If s\ < ei, then one prolongs the integra- 
tion but renormalizes the state vector in the end of the 
first interval at r = : a T + = a T1 / y/\a Tl | 2 + \b Tl | 2 and 

b T + = b Tl j ' \f\a Tl | 2 -I- \b Tl | 2 . If si > ei, then the atom 
emits a spontaneous photon and performs the jump to 
the ground state at r = t\ with A Tl = a T1 = /3 Tl = 0, 
B Tl = 1. Its momentum in the x direction changes for a 
random number from the interval [0, 1] due to the photon 
recoil effect, and the next time step commences. 

We simulate lithium atoms with the relevant transi- 
tion 2S\/2 — 2P3/2, the corresponding wavelength X a = 
670.7 nm, recoil frequency i> lcc = 63 KHz, and the de- 
cay time T sp = 2.73 • 10 -8 s. With the maximal Rabi 
frequency £1o/2tt ~ 126 MHz and the radius of the laser 
beam r = 0.05 cm one gets uj r = 10~ 3 , ay = 400, and 
7 = 0.05. Simulated trajectories in the real space for 50 
spontaneously emitting atoms under the same conditions 
as in Fig.[?]arc shown in Fig.[S] Even though deterministic 
Hamiltonian chaos is masked by random events of spon- 
taneous emission, nevertheless the spatial and momentum 
(not shown) distributions are much broader at those val- 
ues of A at which the Hamiltonian center-of-mass motion 
is chaotic. Namely the chaotic Hamiltonian walking is 
eventually responsible for divergency of atomic beams in 
the real and momentum spaces. 

To simulate a real experiment we consider a beam of 
10 4 lithium atoms with the initial Gaussian distribution 
(the rms o~ x = a p = 2 and the average values, xq = 0, 
and momentum, pg = 10) and compute their distribution 
at a fixed moment of time. In Fig. [7^, we compare the 
atomic position distributions at r = 1000 for the chaotic 
scattering at A = 0.2 (bold curve) and the regular scatter- 
ing at A = 1 (dashed curve) when neglecting spontaneous 
emission. The difference is evident. In the regime of the 
chaotic scattering at A = 0.2 atoms are distributed more 
or less homogeneously over a large distance of 8 wave- 
lengths along the a?-axis whereas they form a few peaks 
in a much more narrow interval under the conditions of 
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the regular scattering at A = 1. Figure [TJd demonstrates 
the distributions of spontaneously emitting atoms at the 
normalized decay rate 7 = 0.05 under the same conditions 
as in Fig. [7jt. The regularly scattered atoms at A = 1 
(dashed curve) form the contrast atomic relief with the bi- 
furcated peaks around the first few SW nodes at x = ±1/4, 
x = ±3/4 and x = 5/4. The distribution of chaotically 
scattered atoms at A = 0.2 (bold curve) has the peaks 
without any bifurcation at x = ±1/4 and x = 3/4 with a 
smaller number of atoms in each one. Moreover, this dis- 
tribution is less contrast as compared to the previous one. 
Thus, we predict that under the conditions of chaotic scat- 
tering there should appear less contrast and more broad- 
ened atomic reliefs as compared to the case of regular scat- 
tering because a large number of atoms are expected to be 
deposited between the nodes as a result of chaotic walk- 
ing along the SW axis. The effect is expected to be more 
prominent under the coherent evolution but it seems to be 
observable with spontaneously emitting atoms as well. 

We predict that experiments on the scattering of atomic 
beams at a SW laser field can directly image chaotic walk- 
ing of atoms along the SW. In a real experiment the final 
spatial distribution can be recorded via fluorescence or ab- 
sorption imaging on a CCD, commonly used methods in 
atom optics experiments yielding information on the num- 
ber of atoms and the cloud's spatial size. The other possi- 
bility is a nanofabrication where the atoms after the inter- 
action with the SW are deposited on a silicon substrate in 
a high vacuum chamber. In this case the spatial distribu- 
tion can be analyzed with an atomic force microscope. As 
to the momentum distribution, it can be measured, for ex- 
ample, by a time-of-flight technique. The modern tools of 
atom optics enable to create narrow initial atomic distri- 
butions in position and momentum, reduce coupling to the 
environment and technical noise, create one-dimensional 
optical potentials, and to measure spatial and momentum 
distributions with high sensitivity and accuracy 



Conclusion. — We have simulated the new type of 
atomic diffraction at a SW without any modulation of its 
parameters and shown that it can be observed in real ex- 
periments. That would be the prove of existence of the 
novel type of atomic motion, chaotic walking in a deter- 
ministic environment. The effect could be used in optical 
nanolithography to fabricate complex atomic structures 
on substrates. 
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